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1 Introduction 



One of the most fundamental challenges pervading theoretical nuclear physics since 
half a century is to understand the properties of nuclei in terms of the underlying 
interactions between the constituents. 

Historically, the first attempt was made by Heisenberg's student H. Euler who 
calculated the properties of nuclear matter in second order perturbation theory 
assuming nucleons interacting via a two-body potential of Gaussian shape. When the 
singular nature of the nuclear potential at short distances ('hard core') was realized, it 
became apparent that conventional perturbation theory is inadequate. Special many- 
body methods had to be worked out. Brueckner and coworkers initiated a method 
which was further developed by Bethe and Goldstone 0]. Alternatively, Jastrow 
suggested to take a variational approach to the nuclear many-body problem. 

In the 1960's, substantial advances in the physical understanding of Brueckner 
theory were made due to the work by Bethe and coworkers (see e. g. the review 
by Day 0). Systematic calculations of the properties of nuclear matter applying 
Brueckner theory started in the late 1960's and continued through the 1970's [7-9] 
(see Ref. |]T^ for a more recent summary). 

The results obtained using a variety of nucleon-nucleon (NN) potentials show a 
systematic behaviour: The predictions for nuclear matter saturation are located along 
a band which does not meet the empirical area, see Fig. 1 and Table 1. (Various 
semi-empirical sources suggest nuclear matter saturation to occur at an energy per 
particle EjA = — 16±1 MeV and a density which is equivalent to a Fermi momentum 
of kp = 1.35 ± 0.05 fm~^ j^, ^.) This phenomenon is denoted by the "Coester band 



)5 



24 . The essential parameter of the Coester band is the strength of the nuclear tensor 
force as measured by the predicted D-state probability of the deuteron or as expressed 
in terms of the wound integral in nuclear matter (see Table 1) [ p!0[| . 

The Brueckner-Goldstone expansion is believed to be convergent in terms of the 



number of hole lines. Calculations by Day |T^, [T^ have confirmed this for the case 



of some realistic potentials. However, three- and four-hole line diagrams contribute 
about 5-7 MeV to the binding energy per nucleon at saturation (cf. Table 1) and, 
thus, are not negligible. In Fig. 1 open symbols represent results obtained in the 
two-hole line approximation, symbols with a cross denote results including the con- 
tributions from three and four hole-lines. It is seen that taking into account up to four 
hole lines leads, indeed, to an improved Coester band as compared to the two-hole 
line approximation; however, the improvement is insufficient to explain the empirical 
saturation point [|^. Results based on the variational approach are in fair agreement 
with Brueckner theory predictions and, thus, also fail to quantitatively explain 
nuclear saturation. 

Since the mid 1970's, there have been comprehensive efforts to check Brueckner 
theory we mention here, in particular, the work using hypernetted chain theory 
Based on this work, there are indications that the two hole-line approxi- 
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Figure 1: Nuclear matter saturation as predicted by a variety of NN potentials (cf. 
Table 1). Open symbols are saturation points obtained in the two- hole line ap- 
proximation, symbols with a cross denote corresponding predictions with three- and 
four-hole lines included, %-numbers refer to the D-state probability of the deuteron 
as predicted by the corresponding potential. The shaded area denotes the empirical 
value. 



mation of Brueckner-Bethe theory may not necessarily be correct. Some published 
variational calculations are in fair agreement with Brueckner-Bethe results if three 
hole-lines plus a four hole-line estimate are included [0]. 

Approaches discussed so far are based on the simplest model for the atomic nu- 
cleus: Nucleons obeying the non-relativistic Schrodinger equation interact through 
a two-body potential that fits low-energy NN scattering data and the properties of 
the deuteron. The failure of this model to explain nuclear saturation indicates that 
we may have to extend the model. One possibility is to include degrees of freedom 
other than the nucleon. The meson theory of the nuclear force suggests to consider, 
particularly, meson and isobar degrees of freedom. Characteristically, these degrees of 
freedom lead to medium effects on the nuclear force when inserted into the many-body 
problem as well as many-nucleon force contributions (see Ref. []TU[ for a comprehen- 
sive review on this subject). In general, the medium effects are repulsive whereas the 
many-nucleon force contributions are attractive. Thus, there are large cancellations 
and the net result is very small. The density dependence of these effects/contributions 
is such that the saturation properties of nuclear matter are not improved [|T0[] . 

We note that the discussion of many-body force effects in the previous paragraph 
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Table 1: Nuclear matter saturation properties as predicted by various NN potentials. 



Potential Ref.'^ Pd (%) k (%) £/A (MeV) kp (fm-i) Ref.^ 



HJ 


11 


7.0 


21 


-7.2 


1.27 


12 


BJ 


13 


6.6 




-8.5 [-14.2] 


1.36 [1.48] 


12 


RSC 


14 


6.5 


14 


-10.3^ [-17.3] 


1.40^= [1.52] 


15 [12] 


Vi4 


16 


6.1 


12 [19] 


-10.8 [-17.8] 


1.47 [1.62] 


17 


Paris 


18 


5.8 


11 


-11.2 [-17.7] 


1.51 [1.63] 


18 [17] 


HMl 


19 


5.8 


11 


-11.8 [-16.9] 


1.48 [1.56] 


19 [17] 


Sch 


20 


4.9 


8 


-20.2 


1.85 


19 


UNG 


21 


4.4 


5 


-23.3 


1.87 


19 


Bonn C 


10 


5.6 


8.1 


-12.1 


1.54 


10 


Bonn B 


10 


5.0 


6.6 


-14.0 


1.61 


10 


Bonn A 


10 


4.4 


5.4 


-17.1 


1.74 
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Given are the saturation energy per nucleon, £/A, and Fermi momentum, 
/cf, as obtained in the two- hole line approximation using the standard 
choice for the single particle potential. Results including three- and four- 
hole line contributions are given in square brakets. The wound integral 
K is given at kp — 1.35 fm^^. Pd is the predicted %-D state of the deuteron. 
" References to the potentials. 
^ References for the nuclear matter calculations. 
Using OBEP for J > 3. 



applies to approaches in which two- and many-body forces are treated on an equal 
footing; i. e. both categories of forces are based on the same meson-baryon interactions 
and are treated consistently. The situation is different if the three-body force is 
introduced on an ad hoc basis with the purpose to fit the empirical nuclear matter 



saturation Such a three-body potential may be large, particularly, if the two- 



body force used substantially underbinds nuclear matter. By construction such a 
three-body force improves the equation of state of nuclear matter as well as the 
description of light nuclei. 

In the 1970's, a relativistic approach to nuclear structure was developed by Miller 
and Green E^]. They studied a Dirac-Hartree model for the groundstate of nuclei 



which was able to reproduce the binding energies, the root-mean-square radii, and the 
single-particle levels, particularly, the spin-orbit splittings. Their potential consisted 
of a strong (attractive) scalar and (repulsive) vector component. The Dirac-Hartree (- 
Fock) model was further developed by Brockmann [^] and by Horowitz and Serot 
pO| , |31| . At about that same time, Clark and coworkers applied a Dirac equation 
containing a scalar and vector field to proton-nucleus scattering [^. The most sig- 
nificant result of this Dirac phenomenology is the quantitative fit of spin observables 
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which were only poorly described by the Schrodinger equation [^. This success and 



Walecka's theory on highly condensed matter ||3J] made relativistic approaches very 
popular. 

Inspired by these achievements, a relativistic extension of Brueckner theory has 



been suggested by Shakin and collaborators |36[, frequently called the Dirac- 
Brueckner approach. The advantage of a Brueckner theory is that the free NN inter- 
action is used; thus, there are no parameters in the force which are adjusted in the 
many-body problem. The essential point of the Dirac-Brueckner approach is to use 
the Dirac equation for the single particle motion in nuclear matter. In the work done 
by the Brooklyn group the relativistic effect is calculated in first order perturbation 
theory. This approximation is avoided and a full self-consistency of the relativistic 
single-particle energies and wave functions is performed in the subsequent work by 



Brockmann and Machleidt |37, 38, 3S], and by ter Haar and Malfliet [40|. Formal 



aspects involved in the derivation of the relativistic G-matrix, have been discussed in 



detail by Horowitz and Serot |41, 42 



The common feature of all Dirac-Brueckner results is that a (repulsive) relativis- 
tic many-body effect is obtained which is strongly density-dependent such that the 
empirical nuclear matter saturation can be explained. In most calculations a one- 
boson-exchange (OBE) potential is used for the nuclear force. In Ref. a more 
realistic approach to the NN interaction is taken applying an explicit model for the 
27r-exchange that involves A-isobars, thus, avoiding the fictitious a boson typically 
used in OBE models to provide intermediate range attraction. It is found in Ref. |^ 



that the relativistic effect for the more realistic model is almost exactly the same as 
that obtained for OBE potentials. This finding seems to justify the use of the OBE 
model. 

It is the purpose of this chapter to present a thorough introduction into the Dirac 
Brueckner approach including the mathematical details of the formalism involved. 
Furthermore, we will present results for nuclear matter, NN scattering in the nuclear 
medium, and finite nuclei. 

2 Relativistic Two-Nucleon Scattering 
2.1 Covariant Equations 

Two-nucleon scattering is described covariantly by the Bethe-Salpeter (BS) equation 
1^. In operator notation it may be written as 



M = V + VgM (1) 

with Ai the invariant amplitude for the two-nucleon scattering process, V the sum 
of all connected two-particle irreducible diagrams and Q the relativistic two-nucleon 



propagator. As this four- dimensional integral equation is very difficult to solve 
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so-called three-dimensional reductions have been proposed, which are more amenable 
to numerical solution. Furthermore, it has been shown by Gross |^6| that the full 



BS equation in ladder approximation does not generate the desired one-body equa- 
tion in the one-body limit (i. e., when one of the particles becomes very massive) 
in contrast to a large family of three-dimensional quasi-potential equations. These 
approximations to the BS equation are also covariant and satisfy relativistic elastic 
unitarity. However, the three-dimensional reduction is not unique, and in principal 
infinitely many choices exist Typically they are derived by replacing Eq. (1) by 
two coupled equations: 

M = W + WgM (2) 

and 

>V = V + V(^-^)W (3) 

where g is a. covariant three-dimensional propagator with the same elastic unitarity 
cut as Q in the physical region. In general, the second term on the r.h.s. of Eq. (3) is 
dropped to arrive at a real simplification of the problem. 

Explicitly, we can write the BS equation for an arbitrary frame (notation and 
conventions of Ref. 



Miq', q\P) = V(g', q\P) + / d^kV{q\ k\P)g{k\P)M{k, q\P) 



with 



g{k\p) 



i2n) 



iP + ky-M^+ie 



2 

(1) 



l(2) 



\P - ky - A/P + ie 



(2) 



(4) 



(5) 



(6) 



where q, k, and q' are the initial, intermediate, and final relative four-momenta, 
respectively (with e. g. A; = {kQ,k)), and P = {Po,P) is the total four-momentum; 

= •y^k^. The superscripts refer to particle (1) and (2). In general, we suppress spin 
(or helicity) and isospin indices. 

Q and g have the same discontinuity across the right hand cut, if 



Im^(A;|P) 



(2vr)4 



2 



(1) 



P + k 



2 



1 



M 



(2) 



P-k 



= lmg{k\P) (7) 

with 5^~^^ indicating that only the positive energy root of the argument of the 5- 
function is to be included. From this follows easily 



Im^(A;|P) 



1 



(1) 



M 



(2) 



5 



5 (iPo + ko- ^ip+k) 5 (i^o -ko- ^^ip-k 



ip+k ^|p-k 



with ^ip±k = + (iP ± k)2. Using the equahty 



5 o^o + fco - ^ip+k )5[-Po-ko- ^ip_k 



5 Po-^i 



P+k 



El 



P k 



5 Uo 



1 



-El 



the imaginary part of the propagator g{k\P) can now be written 



-El 



P k 



(9) 



Im g{k\P) 



1 



M2 



8vr2 i^ip+k ^iP-k 



A 



(1) 



-P + k A 



xS {^Po- -El p+k - -Eip-k 

2 2 



where 



A?(P) 



xS i ko- 



2 

p+k "T o iP-k 



,(2) 



P-k 



2 2' 



/ 7°Ep-7-P + M ' 
^ m(p, Ai)M(p, A,) 



(10) 



(11) 

(12) 



represents the positive-energy projection operator for nucleon i {i = 1, 2) with u{p) 
a positive-energy Dirac spinor of momentum p; Aj denotes either the hehcity or the 
spin projection of the respective nucleon, and Ep = a/M^ + p^. The projection opera- 
tors imply that contributions involving virtual anti-nucleon intermediate states ('pair 
terms') are suppressed. It has been shown in Refs. [49, 50] that these contributions 
are small when the pseudovector coupling is used for the pion. 

Note that lm.g{k\P) is covariant, since Im g{k\P) = Im^(A;|P). 

Using 6{Po - E)= 2E6{s - E^ + P^), where E = Eip+k + ^ip-k and s = P^ = 
Pq — P^, Eq. (10) can be re- written as 



Im5f( 



M2 2 {e. 



|p+k 



+ Ei 



p-k 



^ip+k^ip-k 



A 



(1) 



-P + k A 



P-k 



x5 
xS [ko 



^ip+k + ^ip-k 



-El 

2 2 



p+k 



-El 

2 2 



p-k 



(13) 
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Knowing its imaginary part, we construct g{k\s) by a dispersion integral 

g{k\s) = - r , . l^g{k\s') (14) 

vr Jam^ s' — s — It 

Inserting Eq. (13) in Eq. (14) and integrating, we obtain 

A;^(iP + k)Aj)(iP-k 
(^iP+k + ^ip-k)'-P'-s-^e 
x5( A;o - l^ip+k + i^ip-k ) (15) 
This three-dimensional propagator is known as the Blankenbecler-Sugar (BbS) choice 



5T| , |52[| . By construction, the propagator g has the same discontinuity across the 



right-hand cut as Q; therefore, it preserves the unitarity relation satisfied by Ai. 

Using the angle averages (|P±k)^ ^ |P^-|-k^ and (|P±q)^ ^ |P^ + q^, which 
should be a very good approximation, Eq. (15) assumes the much simpler form 



1 M'^ A« (IP + k) A? (IP - k) 

l^/TTJ i^lp+k ^|P+q ^iP+k + ^^ 

where we used s = 'iE\„ , — P^. 

^P+q 

Assuming W = V, the reduced Bethe-Salpeter equation, Eq. (2), is obtained in 
explicit form by replacing in Eq. (4) Q hy g oi Eq. (16), yielding 

d^k M2 



/fi A 
^V(q',k,.,^^ 



P+k 



Ai^) f iP + k) Af (iP 



X ¥2 '\ ^ -^(k,q|P) (17) 



in which both nucleons in intermediate states are equally far off their mass shell. 

Taking matrix elements between positive-energy spinors yields an equation for the 
scattering amplitude in an arbitrary frame 



d^k , _ M' 



T(q',q|P) = nq',q) + /(^nq',k)^^ 



2 



P+k 



^|p+q - ^|p+k + 



T(k,q|P) 
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where we used 



u 



-i(iP + q')M2(|P - q')V(q', q|P)ni(iP + q)u2(iP - q) 



«i(q')«2(-q')V(q', q)Mi(q)M2(-q) 

^\/(q',q) (19) 

since this is a Lorentz scalar. An analogous statement applies to T. Calculations of 
nuclear matter and of finite nuclei are performed in the rest frames of these systems. 
Thus, Eq. (18) with the necessary medium modifications would be appropriate for 
the evaluation of the nuclear matter reaction (G) matrix. 

In the two-nucleon cm. frame ( i. e. for P = 0), the BbS propagator, Eq. (16), 
reduces to 

„M) ^ W (20) 



le 



which implies the scattering equation 

(2?f''«>''')b:^f3ft«^(''.<i) (21) 

Two-nucleon scattering is considered most conveniently in the two-nucleon cm. frame; 
thus, for calculations of free space two-nucleon scattering in the BbS approximation, 
one would use Eq. (21). 
Defining 



f(q',q) = ^^T(q',q)y'^ (22) 



and 



V(q',q) = ^-^V(q',q)^^, (23) 



which has become known as "minimal relativity" [Q, we can rewrite Eq. (21) as 

^ (q', q) = \>(q', q) + / d'k f (q', k) , f . . f (k, q) (24) 

which has the form of the non-relativistic Lippmann-Schwinger equation. A poten- 
tial defined within an equation which is formally identical to the (non-relativistic) 
Lippmann-Schwinger equation, can then be applied in conventional (non-relativistic) 
nuclear structure physics. This is the practical relevance of Eqs. (22-24). On the 
other hand, for relativistic nuclear structure physics, Eq. (21) is the starting point. 

The BbS propagator is the most widely used approximation. Another choice, 
that has been frequently applied, is the version suggested by Thompson [Q. The 
manifestly covariant form of Thompson's propagator is the same as Eq. (14), but with 
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/4M2 ds' / {s' — s — ie) replaced by d^/s'/ {sfs' — y/s — it). After taking the angle 
average for | |P ± k | and | |P ± q |, this propagator reads 

1 Ai'^ (iP + k) A?) (ip - k) 

The equation for the scattering amplitude in an arbitrary frame is 

d'^k M2 



T(q',q|P) = r(q',q)+ / 7^^(q',k) 



(27r)3 ^Eki?ip+k 

^2E^3i^^(k>q|P) (26) 

For calculations in the rest frame of nuclear matter or finite nuclei, this equation 
together with the necessary medium modifications (M i — > M, Pauli projector Q) 
is appropriate (see Eq. (95) below). In our actual calculations in nuclear matter, 
we replace by -Eip+k and by Eip_^_^ in the denominator of Eq. (26). This 
replacement makes possible an interpretation of the energy denominator in terms 
of differences between single particle energies which are typically defined in the rest 
frame of the many-body system. This allows for a consistent application of this 
equation in nuclear matter and finite nuclei |]^ . The change of the numerical results 
by this replacement is negligibly small (less than 0.1 MeV for the energy per nucleon 
at nuclear matter denisty), since the factor -^^^ is slightly reduced, while the 

term (2i5q — 2ii^k + ^e)~^ is slightly enhanced. 

In the two-nucleon cm. frame (P = 0), the Thompson equation is 

r d^k 1 
T(q'.q) = V(q'.q) + / (Sp^W- '')^2^3I^^(''. (2^) 

This equation is usefull for free-space two-nucleon scattering. In the calculations of 
this chapter, always Thompson's equations are used. 

Again, we may introduce some special definitions. In the case of the Thompson 
equation, it is convenient to define 

M M 
T(q',q) = — r(q',q)- (28) 

-C/g' J^q 

and 

V(q',q) = ^nq',q)|^. (29) 
With this, we can rewrite Eq. (27) as 

T(q', q) = V^(q', q) + / d'k V^(q', k)^^— i^-^f (k, q) (30) 
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which looks hke a Lippmann-Schwinger equation with relativistic energies in the prop- 
agator. 

Since both nucleons are equally off-shell in the BbS or Thompson equation, the 
exchanged bosons tranfer three-momentum only, i. e. the meson propagator is (for a 
scalar exchange) 

(31) 



[q' - q)2 - ml' 



this is also referred to as a static (or non-retarded) propagator. 

Many more choices for g have been suggested in the literature. For most of the 
other choices, the three-dimensional two-nucleon propagator is either of the BbS or 
the Thompson form. However, there may be differences in the 6 function in Eqs. (20) 
and (25). For example, Schierholz |2^ and Erkelenz ||56| propose 6{kQ + ^^/s — E^), 



restricting one particle to its mass shell. It implies the meson propagator 

(32) 



(Eq, -Eq)2-(q'-q)2-m 



2 



However, it has been shown [S^] (see also Appendix E.l of Ref. [^) that the term 
{Eq>-E^y in this propagator has an effect which is opposite to the one obtained when 
treating meson retardation properly. Also, the medium effect on meson propagation in 
nuclear matter which, when calculated correctly, is repulsive |B3 , comes out attractive 



when using the form Eq. (32) Thus, in spite of its suggestive appearance and in 
spite of early believes, the meson propagator Eq. (32) has nothing to do with genuine 
meson retardation and, therefore, should be discarded. This remark applies to the 
scattering of two particles of equal mass. If one particle is much heavier than the 
other one, it may, however, be appropriate to put one particle (namely, the more 
massive one) on its mass shell. 

A thorough discussion of the Bethe-Salpeter equation and/or a systematic study 
of a large family of possible relativistic three-dimensional reductions can be found in 
Refs. [0, ^ Tjon and coworkers have compared results obtained by solving the 
full four-dimensional Bethe-Salpeter equation applying a full set of OBE diagrams 
with those from the BbS and some other three-dimensional equations; for BbS they 



find only small differences as compared to full BS p9| , |50[] . This is also true for the 
Thompson choice, since it differs little from BbS. 

Notice that the quantity T in Eqs. (21) and (27) is invariant, while T of Eq. (24) 
and T of Eq. (30) are not. T is equivalent to the familiar non- relativistic T-matrix. 

The relationship of the invariant T to the ^'-matrix is 

{p'iP2\S\piP2) = {p'iP'2\PiP2) - 2Tii5^^\p\ + P2 - pi - P2)^T(q', q) 

(33) 

with Pi the initial and p'^ {i = 1, 2) the final four momenta of the two interacting 
nucleons. The normalization is {p'\p) = S^^\p' — p). 
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2.2 R-Matrix Equation and Helicity State Basis 

On a computer, real analysis is much faster than complex analysis. It is therefore 
desirable to deal with real quantities whenever possible. The scattering amplitude 
below particle production threshold can be expressed in terms of the real i?-matrix 
(better known as the 'A'-matrix') which is defined by 

T = R- i7rR6{E - Hq)T (34) 

The equation for the real i?-matrix corresponding to the complex T-matrix of Eq. (24) 
is 

^(q, q) = V^(q', q) + ^ / d'k V(q', k) /?(k, q) (35) 

where V denotes the principal value. 

Now, we also need to consider the spin of the nucleons explicitly. The easiest 
way to treat the spin-projections of spin-| particles in a covariant way is to use the 
helicity representation. In our further developments, we will therefore use a helicity 



state basis (cf. Appendix C of Ref. which is based upon Refs. 



The helicity Aj of particle i (with i = 1 or 2) is the eigenvalue of the helicity 
operator ■ Pi/|pi| which is ±|. 

Using hehcity states, the /2-matrix equation reads, after partial wave decomposi- 
tion, 

{\\\'2\R\q\q)\\,\2) = {KK\V\q\q)\\i\2) 

roc _ M 



fOO Wl 

+ dke^—{x[x',\v'iq\k)\hh2) 

huh2 ° q - K 



X {hih2\R^ik,q)\\i\2) (36) 

where J denotes the total angular momentum of the two nucleons. Here and through- 
out the rest of this chapter, momenta denoted by non-bold letters are the magnitudes 
of three-momenta, e. g. q = \q\, k = |k|, etc.; hi and /i2 are the helicities in interme- 
diate states for nucleon 1 and 2, respectively. Equation ( p^ ) is a system of coupled 
integral equations which needs to be solved to obtain the desired matrix elements of 

Ignoring anti-particles, there are 4 x 4 = 16 helicity amplitudes for R'^ . However, 
time-reversal invariance, parity conservation, and the fact that we are dealing with 
two identical fermions imply that only six amplitudes are independent. For these six 
amplitudes, we choose the following set: 



R({q',q) ^ 


{+ + \R'{q\q)\ + 


+) 


Ri{q',q) = 


{+ + \Rj{q',q)\- 


-) 


Riiq',q) = 


{+-\R^{q',q)\ + 


-) 


Ri{q',q) = 


{+-\R^{q',q)\- 


+) 



(37) 
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Riiq'.q) = {+ + \R\q\q)\ + -) 
Ri{q',q) = {+-\R'{q',q)\ + +) 

where ± stands for ±|. Notice that 

Ri{q',q) = Ri{q,q'). (38) 

We have now six coupled equations. To partially decouple this system, it is usefuU 
to introduce the following linear combinations of helicity amplitudes: 



(39) 



'R' 


— dJ dJ 
= ^1-^2 


'R' 


= R3 ~ Ri 




= Ri + i?2 




= -R3 + -R4 






66^J 





We also introduce corresponding definitions for V"^ . (Of course, analogous definitions 
exist for R, R, V, and V.) Using these definitions, Eq. (^) decouples into the follow- 
ing three sub-systems of integral equations: 
Spin singlet 

'R'iq', q) = V^(g', q)+V dk k^-^^ V^(g', k) "^R^k, q) (40) 

JO q^ — /c^ 

Uncoupled spin triplet 

roo M 



/■oo M 

'R\q', q) = ^V\q', q) +V dk k'^— 'V\q' , k) 'R\k, q) (41) 

Jo q'^ — /e^ 

Coupled triplet states 



A/I ^ 

''R\q\q) = ''V\q\q)+V dk (q' , k) R\k, q) 

Jo — /c^ 

+ ^^V^iq',k)^''R\k,q)] 

^^R\q\q) = ^^V\q',q)+V dk e^^—[^^V' (q' , k) ^^R' (k, q) 

Jo q — k, 

+ ^^y-^(g',A;)^^i2^(fc,g)] 

''R\q',q) = ^^V'{q',q)+V dkk^- -['^V' {q' , k) ''''R' {k, q) 

Jo q — 

+ ''V-'{q',k)'^R'{k,q)] 

M ^ ^ 
''R'{q',q) = ''V'iq',q)+V dk e^—['^V' {q' , k) R' {k, q) 

Jo q — /c 

+ ^""V^iq^k) ^^R\k,q)] 



(42) 
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More common in nuclear physics is the representation of two-nucleon states in 
terms of an \LSJM) basis, where S denotes the total spin, L the total orbital angular 
momentum, and J the total angular momentum with projection M. In this basis, we 
will denote the R matrix elements by Rjf ^ = {L' S J M\R\LS J M) . These are obtained 
from the helicity state matrix elements by the following unitary transformation: 
Spin singlet 



pJO 



Uncoupled spin triplet 
Coupled triplet states 



R 



ji 
J J 



pjl 

^J-l,J-l — 


1 


2J+ U 


pjl 

-"-J+l.J+l — 


1 


2J+l[ 


pjl 

-"■J-1,J+1 — 


1 


2J+l[ 


pjl 

-"-J+l.J-l ~ 


1 


2J + 1 [ 



'R' 



(43) 
(44) 



J ^'R' + (J + 1) '^R' + p{J + 1)( ^'i? 
(J + 1) ^^R^ + J 3^^-^ - p{J + 1)( ""R' + 
J(J + 1)( ^^R^ - ^^R^) - J ^^i?^ + ( J + 1) ^^i?-^ 



55 nJ I 66^J-j 
55 dJ I *56^J^ 



J(J + 1)( ^^R-^ - ^^R^) + (J + 1) '""i?^ - J ''i? 



55 p^ 



(45) 



Analogous transformations exist for V , R, R, V, and V. 

Instead of solving the coupled system in the form Eq. (^) , one can also first apply 
the transformation Eq. (|^) to V and R in Eq. ( ^21) yielding a system of four coupled 
integral equations for R in LSJ representation. 



Ri\{q',q) = Vil{q',q) + V dk 

+ Vf4q',k)Ri\ik,q)] 



M 



Vf4q',k) Ri\{k,q) 



Ri\{q',q) = V:!'_{q',q) + V 



dk k^ 



M 



A;2 



+ V:[l{q\k)Ri\{k,q)] 



Ji 



Ri\{q\q) = Vi^_{q\q) + V 



Ji , 



dk k' 



M 



k^ 



[VlL{q'.k) Ri'_{k,q) 



[Vf^{q',k) Ril{k,q) 



+ V^'_iq',k)R!}_ik,q)] 



Ji 



Ri\{q\q) = V:!i{q',q) + V 



dk k^ 



M 



[V:!'_iq',k) Ri\ik,q) 



Ji 



g2 - k"' 



+ V:!l{q',k)Ri\ik,q)] 



Ji 



(46) 



where we used the abbreviations 



pji — pji pji 

^++ — ^J+i,J+i^ ^ — 



pji pji 



pji pji 
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Conventionally, the coupled triplet channels in NN scattering are considered in this 
form. It can also be obtained by decomposing Eq. ( ^51) directly into LSJ states. In a 
non-relativistic consideration, it is the tensor force which couples triplet states with 
L = J±l. 

So far, we have never mentioned the total isospin of the two-nucleon system, T 
(which is either or 1). The reason for this is simply that T is not an independent 
quantum number. Namely, due to the antisymmetry of the two-fermion state, the 
quantum numbers L, S and T have to fulfill the condition 

Thus, for given L and S*, T is fixed. 



On-Shell R-Matrix and Phase Shifts. Phase shifts are a parametrization of the 
unitary S'-matrix which for uncoupled cases is given by 



S.i 



(4J 



Using the above and Eqs. (^3D and (|3^) in partial wave decomposition, one can relate 
the on-energy-shell i?-matrix to the phase shifts as follows: 
Spin singlet 



tan ^5\Ei 



TT 



lab) 



-qM'R'{q,q) 



Uncoupled spin triplet 



tan ^S-\Ei 



IT 



lab 



qM^R\q,q) 



(49) 



(50) 



For the coupled states, a unitary transformation is needed to diagonalize the two-by- 
two coupled i?-matrix. This requires an additional parameter, known as the 'mixing 
parameter' ej. Using the convention introduced by Blatt and Biedenharn [^], the 
eigenphases for the coupled channels are, in terms of the on-shell i?-matrix. 



tan 6-^{Eiab) 



TT 



qM 



Ri- + Ri^ ± 



R' 



cos 2e J 



(51) 



tan2e j{Ei 



2Ri 



lab 



Ri_ - Ri_ 



All /2-matrix elements in these formulae carry the arguments (g, q) where q denotes 
the cm. on-energy-shell momentum that is related to the energy in the laboratory 
system, Eiab-, by 

Eiab = 2q^/M. (52) 
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An alternative convention for the phase parameters has been introduced by Stapp 
et al. known as 'bar' phase shifts. These are related to the Blatt-Biedenharn 



parameters by 

5i + 5i = Si + 5i 

sm{Si-6^) = tan2ej/tan2ej (53) 

sm{6i - 6l) = sin2ej/sin2ej 

We use the 'bar' convention. 

Using the transformation Eq. the phase parameters for the coupled case 

can also be expressed directly in terms of the helicity-state i?-matrix elements; one 
obtains 



tan SUEiab) = -'^qM 



4 



12^-^ _ 34^J _ 4 /j/ J ^ 1) 55^J 



tan2ej(Ez,b) = -2 



(2J+ l)cos2ej 
'(54) 

J(J+1)(12^^- 34j^J) + 55 
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RJ - 34i?J - A\ J{J + 1) 55 i?-/ 



Effective Range Parameters. For low-energy S'-wave scattering, gcot5 can be 
expanded as a function of q 

q J. 1 1 2 

gcot ~ 1 — rq (55) 



tan 5 a 2 

where a is called the scattering length and r the effective range. 

Rewriting this equation for two different (small) on-shell momenta qi and q2 (with 
^lab ~ '^Qi/^ S'lid phase shifts Si, i = 1,2), we can determine the two unknown 
constants a and r, 

A -31 22_ 

^ _ tan Si tan 82 (crf-\ 

^lab - ^lab 

- = —r ^ ELh r (57) 

a 4 tan& ^ ^ 



and 



with z = 1 or 2. 



Using Thompson's Equation. The Thompson equation, Eq. (27), is solved most 
conveniently in 'check' notation, Eqs. ([ 
i?-matrix are obtained from Eqs. (^0D-( 



)- (|30D . The corresponding equations for the 
\) by replacing 



M 



2Eg — 2Ek 



\{Eq + Ek) 

0^ — /c^ 



(5^ 
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and 

^ I — ^ R (59) 

This is easily understood by comparing Eq. (|2^) with Eq. (^Oj). R is defined in analogy 
to Eq. 

The phase shift relation is 

tan '6'{Eiaf,) = -\qE, ^R\q, q) (60) 
and similarly for the other channels. 



3 One-Boson-Exchange Potentials 

3.1 Interaction Lagrangians and OBE Amplitudes 

We use the following Lagrangians for meson-nucleon coupling 

= -i^i,^^Yijd^^^^'^ (61) 

Trips 

Cs = +gsi^^^^'^ (62) 
= -9vi'iy^t^-j^i^^'''i^{d,^^:^-d^^^;^) (63) 

with if) the nucleon and ip"^J^ the meson fields (notation and conventions as in Ref. ||48|| ). 

For isospin 1 mesons, y^*^"^ is to be replaced by r • with r' (Z = 1, 2, 3) the usual 
Pauli matrices, ps, pv, s, and v denote pseudoscalar, pseudovector, scalar, and vector 
coupling/ field, respectively. 

The one-boson-exchange potential (OBEP) is defined as a sum of one-particle- 
exchange amplitudes of certain bosons with given mass and coupling. We use the six 
non-strange bosons with masses below 1 GeV/c^. Thus, 

VoBEP= E K?^"" (64) 

with TT and rj pseudoscalar, a and S scalar, and p and u vector particles. The contri- 
butions from the iso- vector bosons vr, S and p contain a factor Ti • T2- 
The above Lagrangians imply the following OBE amplitudes:^ 

= 7:rV3 K)l'r^iq' - '?)M«(q, Ai)^2(-q', X'^h'^^iq - q),u{-^. A2) 

l^ZTTj nipg 



^Strictly speaking, we give here the potential defined as i times the Feynman amphtude; further- 
more, there is a factor of i for each vertex and propagator; since — 1, we simply ignore these 
factors of i. 
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/[(q'-q)2 + my 

-^^{u{q\ X'.h^q, Ai)S(-q', y,h'u{-q, X2) 
(ZTT) rripg 

+ [{E' ~ E)/{2M)Yu{^, A;)7'^7%(q, Ai)fl(-q', A^)7^7''«(-q, A2) 
+ [{E' - E)/{2M)] [u(q', A'i)7'u(q, Ai)i2(-q', A'2)7^°u(-q, A2) 
+u(q', A'i)7'7°«(q, Ai)lx(-q', A^)7'u(-q, A2)]} 

/[(q'-q)2+m2 ]; (65) 



q'A'iA^lK<^^^|qAiA2) 



(2^) 

(q'A'iA^lFO^^IqAiA 



2 

^2(q',Al)7.(q,Ai)u(-q',A^)7.(-q,A2)/[(q'-q)'+m2]; (66) 



3 



21 



1 - fv - 

{c/„?2(q', A'i)7^w(q, Ai) + r-^u(q', Ai)cr^^i(g' - g)''u(q, Ai)} 



(27r)3 "^^"^^ ' ^ 2M ' 

x{g„S(-q', A'2)7^^(-q, A2) - , ^>'"'^iq' - q)M-^, A2)} 



/[(q'-q 



1^ + m?, 



^ {{gv + fv)u{q',X[)j^u{(i,Xi) 



~ (27r)3 i-^i'" ^ •'"''"^^ ' "1'' 

" 2^"^^'' + + ~ 7M7")]"(q, Ai)} 

x{(g„ + /i,)u(-q',A'2)7^w(-q, A2) 

-^^(-q', A^)!!?' + g)p + {E' E){g^^ - 7^7°)]«(-q, A2)} 

/[(q'-q)2+m^]. (67) 

Working in the two-nucleon cm. frame, the momenta of the two incoming (outgoing) 



nucleons are q and — q (q' and — q')- E = a/ + q^ and E' = \J + q'^. Using the 
BbS or Thompson equation, the four-momentum transfer between the two nucleons 
is (g' — q) = (0, q' — q). The Dirac equation is apphed repeatedly in the evaluations 
of the pf-coupling, the Gordon identity is used in the case of the f-coupling. 
(Note that in Eq. (|67D, second line from the bottom, the term (g' + q)^ carries /i as 
a swtecript to ensure the correct sign of the space component of that term.) The 
propagator for vector bosons is 

, -g^.u + W - q)M' - (1)1^1^1 /gg^ 
— (q' — q)2 — ml 

where we drop the (g' — q)fj,{q' — g)iy-term which vanishes on-shell, anyhow, since 
the nucleon current is conserved. The off-shell effect of this term was examined in 
Ref. and found to be unimportant. 
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The Dirac spinors in helicity representation are given by 



^"(-c.A.) = ;/^( 4)|A.> (70) 



with 

= Xai , IA2) = X-X2 ) (71) 

where x denotes the conventional Pauh spinor. 
We normahze Dirac spinors covariantly, that is 

w(q,A)M(q,A) = l, (72) 

with u = u^'y^. 

At each meson-nucleon vertex, a form factor is apphed which has the analytical 
form 

^<.iw-q)'i- ( a. '''' 

with m„ the mass of the meson involved, Aq the so-called cutoff mass, and an 
exponent. Thus, the OBE amplitudes Eqs. (|65|)-(|67|) are multiplied by JF^. 

In practice it is desirable to have the potential represented in partial waves, since 
scattering phase shifts are defined in such a representation and nuclear structure 
calculations are conventionally performed in an LSJ basis. We will turn to this in 
the next subsection. 

3.2 Partial Wave Decomposition 

The OBE amplitudes are decomposed into partial waves according to 

{KK\V'{q\q)\XiX2) = 27r £'d(cose) <^a„a;-a^ ^(q'A'iA'^lV^lqAiAa) 

(74) 

where 9 is the angle between q and q'; dj^^,{6) are the conventional reduced rotation 
matrices. 

The d:^j^,{9) are expressed in terms of Legendre polynominals, Pj (cos 6*). The 
following types of integrals will occur repeatedly: 



+1 



If - /■■*7rr-^#-:T^JIW-<i) 



1 (q' — q)2 + m' 



2 oL 
a 
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1 (q' — q)^ + ml 
J + 1 J-i (q' — q)"' + 



where t = cos 6'. Notice that these integrals are functions of q', q, rria, Aq,, and ria- 

We state the final expressions for the partial wave OBE amplitudes in terms of 
the combinations of helicity amplitudes defined in Eq. More details concerning 



their derivation are to be found in Appendix E of Ref. ^8 



Pseudo-scalar bosons {t] and tt meson; for tt apply an additional factor of ti • T2): 
Pseudo-vector coupling (pv) 



"Vli = (76) 



pv ^pv y-"- pv J ' pv 

SiyJ _ r-pCl) /-(O) _ f(0) /-(2)^ 

^pv ~ '--'pv \ ^pv -^J ^pv -^J J 

55yJ _ n P(2) 7-(3) 
^pv ^pi> ^ pv 

66T/J _ p(2) r(3) 

I'pj; "-^pil -C^pt, -fj 



with 



2^ 



and 



= E'E-M^ + {E' -Ef{E'E + 3M'^)/{4M 
= -q'q + q'q{E'-Ef/{AM') (77) 

2 

C^P. = ^ A • (78) 
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For the couphng constant, the following relation holds: 



= /. ^ . (T9) 



Scalar coupling (s) (cr and 5 boson; for 5 apply an additional factor of Ti • T2 ): 



QyJ 

s 




+ F^^ 


lyj 
' s 




+ i^i^^ if 


12yJ 
^ s 


= cAK^'^iT 


+ F^^ I? 


UyJ 

s 




+ F^'^ If 


55t/J 
^ s 


= C.F^lf^ 




my J 







(80) 



with 



and 



1 

^' = t MP ^^^^ 



Ff) = -{E'E + M^) 

F^ = q'q (82) 
= M{E' + E). 

Vector bosons (v) {u and p meson; for p apply an additional factor of Ti • T2 ): 
Vector- vector coupling 

= C,,{2E'E-M^)lf 
^Vi = C,,{E'Elf ^ciqlf) 

^^Vi = C„„(2g'g/f + M2/«) (83) 
""yi = C^A^qif^F'Elf) 
^V/, = -C^^MElf 



with 



o2 1 

C„« = — ^ ■ (84) 



Tensor-tensor coupling 

= c«{(g'' + g')(3£;'£; + M2)7f 
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Table 2: Parameters of the relativistic OBEPs used in this work. 



Potential A Potential B Potential C 

Meson Parameters 

m„(MeV) gl/iT: Ac(GeV) gl/i^n Ac(GeV) gl/i^n A„(GeV) 



TT 


138.03 


14.9 


1.05 


14.6 


1.2 


14.6 


1.3 


»? 


548.8 


7 


1.5 


5 


1.5 


3 


1.5 


P 


769 


0.99 


1.3 


0.95 


1.3 


0.95 


1.3 


U! 


782.6 


20 


1.5 


20 


1.5 


20 


1.5 


s 


983 


0.7709 


2.0 


3.1155 


1.5 


5.0742 


1.5 


c 


550 


8.3141 


2.0 


8.0769 


2.0 


8.0279 


1.8 



Deuteron 
-ea (MeV) 
Pd (%) 
Qd (fm2) 
(mat) 
As (fm-i/2) 

D/S 



2.22459 
4.47 
0.274 " 
0.8543 " 
0.8984 
0.0255 



2.22468 
5.10 
0.279 " 
0.8507 " 
0.8968 
0.0257 



2.22450 
5.53 
0.283 " 
0.8482 " 
0.8971 
0.0260 



[2.224575(9)] 

H 

[0.2860(15)] 
0.857406(1)] 
[0.8846(8)] 
[0.0264(12)] 



Low Energy Scattering 

a„p (fm) -23.752 

rnp (fm) 2.69 

at (fm) 5.482 

n (fm) 1.829 



-23.747 
2.67 
5.474 
1.819 



-23.740 
2.66 
5.475 
1.821 



[-23.748(10)] 
[2.75(5)1 
[5.419(7)] 
[1.754(8)] 



The meson parameters define the potentials. The deuteron and low energy scattering 
parameters are predictions by these potentials. The experimental values are given in square 
brakets. For notation and references to the empirical values see Table 4.2 of Ref. [10].. 
It is always used: fp/gp = 6.1 and fuj/gu = 0.0. 
" Meson-exchange current contributions not included. 



with 



+[q'^ + g2 - 2{3E'E + M'')]q'q I? - 2q'Y if} 
Xi = CuWV + iq" + q'){E'E-M')]lf 
+2{E'E + M^)q'qlf 

-{q'^ + q^ + AE'E)q'q if - 2fq^ if} (85) 
'Xi = Cu{m'-3{q" + q')]q'qlf 

+ [6g'2g2 _ (^'2 ^ g2^^E'E + 3M2)] if + 2{E'E + M^)q'q if} 
'%i = Cu{-{q" + q' + 4E'E)q'qlf-2q'Ylf 

+ [Aq'\^ + (g" + q^E'E - M^)] if + 2{E'E + M^)q'q if} 
= CuM{[E\q" + q') + E{3q"-q')]lf-2iE' + E)q'qlf} 
= CuM{[E{q'' + q') + E\3q'-q")]lf-2{E' + E)q'qlf} 



f 1 

^" ^ 4^ MP ^^^^ 
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Figure 2: Neutron- proton phase shifts as predicted by the three potentials defined 
in Table 2. The solid, dashed, and dotted lines refer to Potentials C, B, and A, 
respectively. The open circles represent the phase shift analysis by Arndt et al. [70]; 
the full dots and the full square are from two analyses by Bugg, Refs. [71] and [72], 
respectively 



Vector-tensor coupling 



with 





= Cyt 








M[-{q''+q^)lf+2q'qlf 




— Cvt 


M[6q'qlP -3{q" + q')I^^] 






M[2q'qlP-{q'' + q')lf^] 


55t/J 




{E'q^ + SEq''') /f 




— Cvt 


{Eq'^ + 3E'q^)lf^ 



/-y _ 9vfv 



(87) 



AnM 27rM2 ' ' 

We use units such that h = c = 1. Energies, masses and momenta are in MeV. 
The potential is in units of MeV^^. The conversion factor is he = 197.327 MeV fm. 
If the user wants to relate our units and conventions to those used in the work of 
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other researchers, he/she should compare our Eq. ( ^01 ) and our phase shift relation 
Eq. (|49|) with the corresponding equations in other works. A computer code for the 
relativistic OBE potential is published in Ref. ||69[| . 

We note that when, in the nuclear medium (nuclear matter), the free Dirac spinors, 
Eq. (69), are replaced by the in- medium Dirac spinors, Eq. (92), then all E and E' 
in the above expressions are replaced by E and E', respectively, and all M by M; an 
exception from this rule are the factors Cu, Eq. (86), and C„t, Eq. (88), which are in 
the medium: ^ 

^ fv ^ 9vfv 1 

" ~ 47rM2 87rM2 ' ~ AnM 2nM'^ ' 

notice that the factor f^/M comes from the Lagrangian, Eq. (63), and we assume 
here that it does not change in the medium. 



3.3 Meson Parameters and Two-Nucleon Predictions 

We give here three examples for relativistic momentum-space OBEPs, which have 
become known as the Bonn A, B, and C potentials[[ They are presented in Appendix 
A, Table A. 2, of Ref. [|l^]. These potentials use the pseudovector coupling for vr and 
7] and are constructed in the framework of the Thompson equation, Eqs. (27-30). 

In Table 2 we give the meson parameters as well as the predictions for the deuteron 
and the low-energy scattering parameters. Phase shifts for neutron-proton scattering 
are displayed in Fig. 2. 



4 The Dirac-Brueckner Approach 

As mentioned in the Introduction, the essential point of the Dirac-Brueckner approach 
is to use the Dirac equation for the single-particle motion in nuclear matter 

(^- Af- [/)m(p,s) = (89) 

or in Hamiltonian form 



(a ■ p [3M + pU)u{p, s) = epu{p, s) (90) 

with 

U = Us + -f''Uv (91) 

where Us is an attractive scalar and Uy (the time-like component of) a repulsive 
vector field. (Notation as in Ref. [48] ; j3 = 7°, a' = 7^7'.) M is the mass of the free 
nucleon. 



^It is customary to denote the OBE parametrizations of the Bonn fuU model by the first 
letters of the alphabet, A, B or C, see e. g. Appendix A of Ref. |10[. 
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The fields, Us and Uy, are in the order of several hundred MeV and strongly 
density dependent (numbers will be given below). In nuclear matter they can be 
determined self-consistently. The resulting fields are in close agreement with those 
obtained in the Dirac phenomenology of scattering. 

The solution of Eq. (89) is 



E„ + M f I \ 



E„+M 



M = M + Us, (93) 



u{p, s) 

with 



Ep = VM2 + p2, (94) 

and Xs a Pauli spinor. The covariant normalization is u{p, s)u{p, s) = 1. Notice that 
the Dirac spinor Eq. (92) is obtained from the free Dirac spinor by replacing M by 
M. 

As in conventional Brueckner theory, the basic quantity in the Dirac-Brueckner 
approach is a G-matrix, which satisfies an integral equation. In this relativistic ap- 
proach, a relativistic three-dimensional equation is chosen. Following the basic phi- 
losophy of traditional Brueckner theory, this equation is applied to nuclear matter in 
strict analogy to free scattering. 

Including the necessary medium effects, the Thompson equation in the nuclear 
matter rest frame reads [cf. Eq. (26)] 



G(q',q|P,5) = V^(q',q)+ / ^^V{^,k) f ^^^f^ g(k,q|P,i) (95) 

l^^J ^ ip+k ^ - ^^iP+k 

with 

~z = 2Eip+^. (96) 

P is the cm. momentum, and q, k and q' are the initial, intermediate and final 
relative momenta, respectively, of the two particles interacting in nuclear matter. 
The Pauli operator Q projects onto unoccupied states. In Eq. (95) we suppressed 
the kp dependence as well as spin (helicity) and isospin indices. For ||P±q| and 
||P ± k| the angle average is used. 

The energy per nucleon in nuclear matter, which is the objective of these calcu- 
lations, is considered in the nuclear matter rest frame. Thus, the G-matrix is needed 
for the nuclear matter rest frame. Equation (95) gives this nuclear matter G-matrix 
directly in that rest frame. Alternatively, one can calculate the G-matrix first in 
the two-nucleon center-of-mass (cm.) system (as customary in calculations of the 
T- matrix in two-nucleon scattering), and then perform a Lorentz transformation to 
the rest frame. This method, which is described in detail in Ref. [Q, is, however. 



complicated, involved, and cumbersome. The advantage of our procedure is that it 
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Figure 3: The repulsive relativistic effect in nuclear matter as obtained in a Dirac- 
Brueckner-Hartree-Fock calculation using Potential B. Saturation points from conven- 
tional calculations are displayed in the background (cf. Fig. 1). The checked rectangle 
represents the empirical value for nuclear matter saturation. 

avoids this complication. Further treatments of Eq. (95) can follow the lines estab- 
lished from conventional Brueckner theory, as e. g. the use of the angle averaged 
Pauli projector etc.. Numerically the equation can be solved by standard methodes 
of momentum space Brueckner calculations 0. 

The essential difference to standard Brueckner theory is the use of the potential 
V in Eq. (95). Indicated by the tilde, this meson-theoretic potential is evaluated 
by using the spinors of Eq. (92) instead of the free spinors applied in scattering as 
well as in conventional ('non-relativistic') Brueckner theory. Since Us (and M) are 
strongly density dependent, so is the potential V. M decreases with density. The 
essential effect in nuclear matter is a suppression of the (attractive) a-exchange; this 
suppression increases with density, providing additional saturation. It turns out (see 
figures below) that this effect is so strongly density- dependent that the empirical 
saturation and incompressibility can be reproduced. Furthermore, the prediction for 
the Landau parameter /o is considerably improved without deteriorating the other 
parameters (see Table 6 below). Note, that sum rules require /o > —1 at nuclear 
matter density [l73| . 

The single-particle potential 



U{m) 



M 



{m\U\m) 




(97) 



'm 
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Table 3: Results of a relativistic Dirac-Brueckner calculation in comparison to the 
corresponding non- relativistic one using Potential B. 



relativistic non-relativistic 



kp £/A M/M Us Uv K £/A M/M Uq" k 

(fm-i) (MeV) (MeV) (MeV) (%) (MeV) (MeV) (%) 



0.8 


-7, 


.02 


0.855 


-136, 


,2 


104, 


.0 


23, 


,1 


-7.40 


0.876 


-33, 


,0 


26 


.5 


0.9 


-8, 


.58 


0.814 


-174, 


,2 


134, 


.1 


18, 


,8 


-9.02 


0.836 


-41, 


.0 


21, 


.6 


1.0 


-10 


.06 


0.774 


-212, 


,2 


164, 


.2 


16, 


,1 


-10.49 


0.797 


-49, 


,0 


18, 


.5 


1.1 


-11, 


.18 


0.732 


-251, 


,3 


195, 


.5 


12, 


,7 


-11.69 


0.760 


-58, 


,1 


14, 


.2 


1.2 


-12, 


.35 


0.691 


-290, 


,4 


225, 


.8 


11, 


,9 


-13.21 


0.725 


-68, 


,5 


12, 


.9 


1.3 


-13 


.35 


0.646 


-332, 


,7 


259, 


.3 


12, 


,5 


-14.91 


0.687 


-80, 


,5 


13 


.1 


1.35 


-13 


.55 


0.621 


-355, 


,9 


278, 


.4 


13, 


,0 


-15.58 


0.664 


-86, 


,8 


13 


.2 


1.4 


-13, 


.53 


0.601 


-374, 


,3 


293 


.4 


13, 


,8 


-16.43 


0.651 


-93, 


,2 


13 


.5 


1.5 


-12, 


.15 


0.559 


-413, 


,6 


328, 


.4 


14, 


,4 


-17.61 


0.618 


-106, 


,1 


13 


.0 


1.6 


-8, 


.46 


0.515 


-455, 


,2 


371, 


.0 


15, 


,8 


-18.14 


0.579 


-119, 


,4 


12, 


.7 


1.7 


-1, 


.61 


0.477 


-491, 


,5 


415, 


.1 


18, 


,4 


-18.25 


0.545 


-133, 


,2 


13, 


.2 


1.8 


+9, 


,42 


0.443 


-523, 


,4 


463 


.6 


21, 


,9 


-17.65 


0.489 


-147, 


,2 


14, 


.3 


1.9 


25, 


,26 


0.418 


-546, 


,7 


513, 


.5 


25, 


,2 


-16.41 


0.480 


-160, 


,7 


15, 


.0 


2.0 


47, 


.56 


0.400 


-563, 


.6 


568, 


.6 


27, 


,5 


-13.82 


0.449 


-173, 


.6 


15, 


.3 


2.1 


77, 


,40 


0.381 


-581, 


.3 


640, 


.9 


30, 


,2 


-9.70 


0.411 


-186, 


.3 


15, 


.7 


2.2 


114, 


.28 


0.370 


-591, 


,2 


723, 


.5 


33, 


,3 


-3.82 


0.373 


-198, 


,1 


16, 


.3 



As a funtion of the Fermi momentum kp, we list: the energy per nucleon £/A, 
M/M, the single-particle scalar and vector potentials Us and Uy, and the wound 
integral k. 

Uo is to be compared to Us + Uy, cf. Eq. (104). 

is the self-energy of the nucleon which is defined in terms of the G-matrix formally 
in the usual way 

U{m) = Re ^ — — —{'mn\G{z)\mn — nm) (98) 

n<kF ^nErn 

where m denotes a state below or above the Fermi surface (continuous choice). The 
constants Us and Uy are determined from Eqs. (97) and (98). Note that Eq. (91) is an 
approximation, since the scalar and vector fields are in general momentum dependent; 
however, it has been shown that this momentum dependence is very weak p8[] . 
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Table 4: Partial wave contributions to the energy in nuclear matter (in MeV) for a 
non-relativistic and a relativistic calculation using Potential B. 



1.1 fm-i 



kp = 1.35 fm 



1.6 fm- 



State 


non-rel. 


relativ. 


non-rel. 


relativ. 


non-rel. 


relativ. 




-10.79 


11.18 


-16.01 


-16.42 


-21.51 


-20.36 




-2.07 


-1.48 


-3.74 


-1.34 


-5.61 


4-2.17 


'Pi 


1.73 


1.77 


3.25 


3.45 


5.33 


6.08 


i 


4.71 


5.27 


9.77 


12.33 


17.69 


26.65 




-15.41 


-14.16 


-20.10 


-17.10 


-23.77 


-17.03 




0.59 


0.57 


1.38 


1.29 


2.64 


2.25 




-0.95 


0.91 


-2.28 


-2.01 


-4.57 


3.39 




-1.70 


1.62 


-4.00 


-3.56 


-7.71 


5.99 


'P2 


-3.10 


-2.92 


-7.06 


-6.28 


-13.31 


-10.73 


'F, 


-0.19 


-0.18 


-0.54 


-0.44 


-1.19 


-0.67 


'F^ 


0.32 


0.31 


0.80 


0.75 


1.60 


1.40 


'Fs 


0.56 


0.55 


1.51 


1.43 


3.20 


2.87 


'Ds 


-0.01 


0.00 


-0.03 


0.00 


-0.11 


-0.02 


'G, 


0.06 


0.06 


0.20 


0.18 


0.49 


0.41 


J > 4 


-0.34 


-0.33 


-1.07 


-0.98 


-2.57 


-2.13 


Total 














potential energy 


-26.61 


-24.25 


-37.93 


-28.72 


-49.38 


-18.51 


Kinetic energy 


14.91 


13.07 


22.35 


15.16 


31.23 


10.05 


Total energy 


-11.69 


-11.18 


-15.58 


-13.55 


-18.14 


-8.46 



The energy per nucleon in nuclear matter is 



]4 = T E -^HTPm+M\m) + — ^ 



In Eqs. (98) and (99) we use 



2^ m^kp EmEn 



Em + E„ 



{mn\G{z)\mn-nm)-M . (99) 



(100) 



The expression for the energy, Eq. (99), is denoted by the Dirac-Brueckner-Hartree- 
Fock (DBHF) approximation. If M is replaced by M, we will speak of the Brueckner- 
Hartree-Fock (BHF) approximation, since this case, qualitatively and quantitatively, 
corresponds to conventional non-relativistic Brueckner theory. Thus, we will occa- 
sionally denote the DBHF calculation by 'relativistic' and the BHF calculation by 
'non-relativistic' (though, strictly speaking, all our calculations are relativistic). 
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Table 5: Energy per particle, £ / A, Fermi momentum, kp, and kompression modulus, 
K, at saturation for nuclear matter with and without relativistic effects. 



relativistic non-relativistic 



Potential 


£/A 
(MeV) 


kp 
(fm-i) 


K 

(MeV) 


£/A 

(MeV) 


kp 

(fm-i) 


K 

(MeV) 


A 


-15.59 


1.40 


290 


-23.55 


1.85 


204 


B 


-13.60 


1.37 


249 


-18.30 


1.66 


160 


C 


-12.26 


1.32 


185 


-15.75 


1.54 


143 



In Equations. (97-99) the states \m) and \n) are represented by Dirac spinors of 
the kind Eq. (92) and an appropriate isospin wavefunction, (m| and (n| are the adjoint 



Dirac spinors u = m^7° with uu = 1; Em = y M"^ + p^. The states of the nucleons 
in nuclear matter, w, are to be normalized by w'^w = 1. This is achieved by defining 



w = yM/E X u which explains factors of M/E in Eqs. (97-99). 

The first term on the r.h.s. of Eq. (99) — the 'kinetic energy' — is in more explicit 
form 

\_ ^ MM + vl ^ ^^^^^ 

m<kp 

The single particle energy is 

M 

em = -^{m\j ■pm + M\m) + U{m) (102) 

Em 

= Em + Uv (103) 
= Em-M + M + Uo. (104) 

5 Results for Nuclear Matter 

We apply now the three OBE potentials presented in Sect. 3.3 to nuclear matter. 
We stress again, that our potentials use for vrA^A^ the coupling of derivative (pseu- 
dovector) type. This is improtant because the pseudoscalar vrA^A^ coupling leads to 
an unrealistcally large attractive medium effect |l38| . 
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Table 6: Landau parameters at various densities from a non-relativistic and a rela- 
tivistic nuclear matter calculation using Potential B. 



(fm-i) 


density 




/o 




30 


9', 


1.0 


0.4po 


relativistic 


-1.37 


0.57 


0.22 


0.66 






non-reiativisic 


-1.50 


0.62 


0.16 


0.66 


1.35 


Po 


relativistic 


-0.79 


0.35 


0.28 


0.67 






non-relativisic 


-1.27 


0.38 


0.15 


0.67 


1.7 


2po 


relativistic 


0.56 


0.29 


0.36 


0.68 






non-relativisic 


-0.99 


0.20 


0.14 


0.69 


2.0 


3.25po 


relativistic 


2.21 


0.37 


0.38 


0.69 






non-relativisic 


-0.71 


0.09 


0.11 


0.71 




(empiri 


ical) 


(0±0.2) 


0.8) 


0.2) 


0.9) 



The main difference between the three potentials applied here, is the strength of 
the tensor force as reflected in the predicted D-state probability of the deuteron, Pq. 
With Pd = 4.5%, Potential A has the weakest tensor force. Potential B and C predict 
5.1% and 5.5%, respectively (cf. Table 2). It is well-known that the strength of 
the tensor force determines the location of the nuclear matter saturation point on the 
Coester band ||2^. To find out the structure of the Coester band, predictions from 
more than one potential are needed. 

All results presented in this section are obtained either in the Bruckner-Hartree- 
Fock (BHF) or the Dirac-Brueckner-Hartree-Fock (DBHF) approximation; as men- 
tioned before, occasionally, we will denote these two methods also as the 'non- 
relativistic' and the 'relativistic' calculation, respectively. 

The repulsive relativistic effect in nuclear matter as created by the DBHF method 
in shown in Fig. 3. In addition, in Table 3, we list the following quantities as a 
function of the Fermi momentum kp'- the energy per nucleon, S/A, the ratio M/M, 
the single-particle scalar and vector potentials Us and Uy, and the wound integral k. 
(For the definition of k and for explicit formulae appropriate for the momentum-space 
framework, see section 5 of Ref. 0.) k can be understood as the expansion parameter 
for the hole-line series. 

As mentioned in Sect. 4, the suppression of the a contribution can be understood 
in simple terms by considering the covariant one-a-exchange amplitude, for q' = 
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Table 7: Nuclear matter results of a relativistic Dirac-Brueckner calculation applying 
Potential A. Notation as in Table 3. 



relativistic non- relativistic 





£/A 


M/M 


Us 


Uv 


K 


£/A 


fm-i) 


(MeV) 




(MeV) 


(MeV) 


(%) 


(MeV) 


0.8 


-7.27 


0.857 


-134.3 


101.6 


22.6 


-7.60 


0.9 


-8.97 


0.817 


-172.1 


131.1 


18.1 


-9.38 


1.0 


-10.62 


0.777 


-209.8 


160.6 


15.2 


-11.01 


1.1 


-11.96 


0.736 


-248.5 


191.1 


11.5 


-12.44 


1.2 


-13.44 


0.692 


-288.8 


222.0 


10.6 


-14.24 


1.3 


-14.86 


0.647 


-331.6 


255.0 


11.0 


-16.35 


1.35 


-15.32 


0.621 


-355.7 


274.7 


11.5 


-17.28 


1.4 


-15.59 


0.601 


-374.9 


289.8 


12.1 


-18.41 


1.5 


-14.88 


0.557 


-416.3 


325.7 


12.7 


-20.25 


1.6 


-11.96 


0.511 


-459.6 


368.6 


14.2 


-21.64 


1.7 


-5.88 


0.470 


-497.2 


412.8 


16.9 


-22.76 


1.8 


+4.44 


0.435 


-530.4 


461.6 


20.4 


-23.38 


1.9 


19.72 


0.409 


-554.8 


512.0 


23.7 


-23.54 


2.0 


41.62 


0.390 


-572.4 


567.5 


26.3 


-22.60 


2.1 


71.20 


0.371 


-590.2 


640.3 


29.1 


-20.42 


2.2 


106.52 


0.366 


-595.0 


732.5 


31.9 


-16.74 



q and Aj = (as used in the Hartree approximation), in which case, due to the 
covariant normalization of the Dirac spinors, the numerator becomes 1. Since the 
physical states of the nucleons in nuclear matter, w, are to be normalized by w'^w = 1 



implying w = yM/E x u, the sigma (as any other) contribution gets the (scalar 
density) factor (M/E)'^ (see second term on the r.h.s. of Eq. (99)) which decreases with 
decreasing M (i. e. increasing density). A corresponding consideration for the time- 
like (70) component of w-exchange would lead to no changes for that contribution. 
However, due to the exchange term and correlations there is a small enhancement of 
the repulsion created by the u with density. 

/^From the numbers given in Table 3 it is seen that the relativistic effect on the 
energy per nucleon, A{S/A)rei (i- e. the difference between the relativistic and non- 
relativistic calculation), is well fitted by the ansatz 

A{S/A)rei ^ 2 MeV X {p/pof\ (105) 



which is suggested by an estimate by Keister and Wiringa 
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Table 8: Nuclear matter results of a relativistic Dirac-Brueckner calculation applying 
Potential C. Notation as in Table 3. 



relativistic non- relativistic 





£/A 


M/M 


Us 


Uv 


K 


£/A 


(fm-i) 


(MeV) 




(MeV) 


(MeV) 


(%) 


(MeV) 


0.8 


-6.80 


0.851 


-140.3 


108.6 


23.0 


-7.21 


0.9 


-8.30 


0.812 


-176.6 


137.2 


19.1 


-8.79 


1.0 


-9.69 


0.773 


-212.9 


165.8 


16.9 


-10.20 


1.1 


-10.64 


0.732 


-252.1 


197.5 


13.6 


-11.15 


1.2 


-11.57 


0.688 


-292.8 


229.8 


13.0 


-12.46 


1.3 


-12.25 


0.644 


-334.0 


262.8 


13.8 


-13.87 


1.35 


-12.24 


0.620 


-356.6 


281.8 


14.5 


-14.35 


1.4 


-11.99 


0.601 


-374.2 


296.4 


15.3 


-15.02 


1.5 


-10.06 


0.561 


-411.8 


330.9 


16.0 


-15.73 


1.6 


-5.72 


0.519 


-451.7 


373.2 


17.5 


-15.65 


1.7 


+ 1.81 


0.482 


-486.3 


416.6 


20.2 


-15.01 


1.8 


13.50 


0.450 


-516.8 


464.6 


23.7 


-13.50 


1.9 


29.85 


0.426 


-538.8 


513.8 


26.9 


-11.22 


2.0 


52.48 


0.409 


-554.6 


567.7 


29.2 


-7.39 


2.1 


82.54 


0.391 


-571.8 


639.4 


31.8 


-1.75 


2.2 


119.89 


0.379 


-583.0 


717.7 


34.9 


+5.88 



Furthermore, Table 3 shows that up to nuclear matter density the wound integral 
n is slightly smaller for the relativistic approach than for the non-relativistic one. This 
implies that in this region the relativistic many-body scheme should be slightly better 
convergent. Beyond nuclear matter density, the situation is reversed. In addition, it 
is amusing to note that for all values of kp the ratio M/M is almost the same for 
the non-relativistic and the relativistic approach. Low values for M/M have often 
been critizised. However, they are not a consequence of the relativistic approach but 
are due to the Brueckner-Hartree-Fock approximation. Higher order corrections will 
enhance this quantity. For a recent discussion of the effective mass in relativistic and 



non-relativistic approaches see Refs. |^ 76 



The representation of nucleons by Dirac spinors with an effective mass, M, can 
be interpreted, as taking virtual nucleon-antinucleon excitations in the many-body 
environment (many-body Z-graphs) effectively into account |7^. This can be made 
plausible by expanding the spinors Eq. (92) in terms of (a complete set of) spinor 
solutions of the free Dirac equation which will necessarily also include solutions rep- 
resenting negative energy (antiparticle) states |^ . 
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kp ( fm ' ) 

Figure 4: Results from calculations with a family of relativistic potentials revealing 
a new Co ester band which meets the empirical area; solid lines: relativistic, dashed 
lines: non-relativistic calculations. For comparison, saturation points from conven- 
tional calculations are displayed in the background (cf. Fig. 1). The shaded square 
denotes the empirical value for nuclear matter saturation. 

In Table 4, we compare the contributions in various partial- wave states as ob- 
tained in a relativistic calculation to that from the corresponding non-relativistic 
one. Detailed investigations have shown that the repulsive relativistic effect seen in 
that table for the P-wave contributions is essentially due to a suppression together 
with a signature of spin-orbit force enhancement. The change of the ^5*0 contribution 
is so small, because of a cancelation of effects due to cr and p. Apart from a reduc- 
tion, the repulsive effect in ^Si is due to a suppression of the twice iterated one-pion 
exchange for reasons quite analogous to the sigma suppression. 

A comparison between relativistic and non-relativistic Brueckner-Hartree-Fock 
calculations for all three potentials is shown in Fig. 4. For the non- relativistic ap- 
proach, the three saturation points are clearly on the Coester band. However, using 
the relativistic method, the saturation points are located on a new band which is 
shifted towards lower Fermi momenta (densities) and even meets the empirical area. 
This is a very desirable effect. The reason for this shift of the Coester band is the 
additional strongly density dependent repulsion which the relativistic approach gives 
rise to. In Table 5, the saturation points for the different potentials are given for 
the relativistic as well as the non-relativistic calculation. In the relativistic case, the 
incompressibihty of nuclear matter using Potential B is about 250 MeV which is in 
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Figure 5: In-medium np differential cross sections at 50 MeV laboratory kinetic en- 
ergy, as obtained for various densities. 



satisfactory agreement with the empirical value of 210 ± 30 MeV |7^. Note that in 
the relativistic Walecka model, 540 MeV is obtained for the compression modulus 

For completeness, we present in Table 6 the Landau parameters at various densities 
as obtained in a relativistic as well as a non-relativistic nuclear matter calculation. 
Based on the nuclear matter G-matrix, the effective particle-hole interaction at the 
Fermi surface is calculated and, multiplied by the density of states kpM / {h^T{'^), 
parametrized by: F = f + f'ri ■ T2 + gcri ■ + g'ox ■ (^2^1 ■ '^2- /,From an expansion 
of the parameters in terms of Legendre polynominals. Pi, we give in Table 6 the 
coefficient for / = 0. For more details and the empirical values see Ref. Notice 
that the prediction for the Landau parameter /o is considerably improved in the 
relativistic calculation without deteriorating the other parameters. Sum rules require 
/o > —1 at nuclear matter density Finally, Tables 7 and 8 contain more results 



for Potentials A and C. 

Concerning nuclear matter at higher densities and neutron matter, we like to refer 
the interested reader to Ref. [l76|. 



6 In-Medium NN Cross Sections 

The investigation of in-medium NN scattering is of interest for intermediate-energy 
heavy-ion reactions. Experimentally, nucleus-nucleus collisions at intermediate ener- 
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Figure 6: In-medium np differential cross sections at 300 MeV laboratory k;inetic en- 
ergy, as obtained for various densities using the Bonn A potential (solid, long-dashed, 
and short-dashed lines). The dotted line shows the parametrisation by Cugnon et al. 
[92]. 



gies provide a unique opportunity to form a piece of nuclear matter in the laboratory 
with a density up to 2-3po (with po? in the range of 0.15 to 0.19 fm~^, the saturation 
density of normal nuclear matter; in this section we use po=0.18 fm~^) Thus 



it is possible to study the properties of hadrons in a dense medium. Since this piece 
of dense nuclear matter exists only for a very short time (typically 10~^^ -10~^^ s), 
it is necessary to use transport models to simulate the entire collision process and 
to deduce the properties of the intermediate stage from the known initial conditions 
and the final-state observables. At intermediate energies, both the mean field and 
the two-body collisions play an equally important role in the dynamical evolution 
of the colliding system; they have to be taken into account in the transport models 
on an equal footing, together with a proper treatment of the Pauli blocking for the 
in-medium two-body collisions. The Boltzmann-Uehling-Uhlenbeck (BUU) equation 
5^ , Q and quantum molecular dynamics (QMD) |5S[, as well as their relativistic 



extensions (RBUU and RQMD) |^6|, |8^, |8^, are promising transport models for 
the description of intermediate-energy heavy-ion reactions. In addition to the mean 
field [Q, the in- medium NN cross sections are also important ingredients of these 



transport models. Specifically, in-medium total as well as differential NN cross sec- 
tions are needed by these models in dealing with the in-medium NN scattering using 
a Monte Carlo method. 

In this section, we will calculate the elastic in-medium NN cross sections in a 
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Figure 7: In-medium np total cross sections as obtained for various densities using 
the Bonn A potential (solid, long-dashed, and short-dashed lines). The dotted line 
shows the parametrisation by Cugnon et al. [92]. 



microscopic way. We base our investigation on the Bonn potentials presented in 
Sect. 3 and the Dirac-Brueckner approach explained in Sect. 4. The original, detailed 



paper on this subject is Ref. |]91 



In-medium NN cross sections can be calculated from the in-medium G matrix. 
Here, we evaluate the G-matrix in the center-of-mass (cm.) frame of the two inter- 
acting nucleons, i. e., we use Eq. (95) with P = 0. For the starting energy, Eq. (96), 
we have now: z = 2Eg = 2y/m^ + q'^, where q is related to the kinetic energy of the 
incident nucleon in the "laboratory system" (Eiab), in which the other nucleon is at 
rest, by: Eiat = 2q^ /m. Thus, we consider two colliding nucleons in nuclear matter. 
The Pauli projector is represented by one Fermi sphere as in conventional nuclear 
matter calculations. This Pauli-projector, which is originally defined in the nuclear 
matter rest frame, must be boosted to the cm. frame of the two interacting nucleons. 



For a detailed discussion of this and the explicit formulae, see Refs. |i2 |. 

In Fig. 5 we show the results for the in-medium np differential cross sections as 
a function of the cm. angle at Eiab=50 MeV and in Fig. 6 at 300 MeV. We consider 
the medium densities p=0 (free-space scattering), (l/2)po, Po, and 2po- The results 
are obtained by using the Bonn A potential. At low incident energies (Fig. 5), the 
np differential cross section always decreases with increasing density, at both forward 
and backward angles. At high incident energies (Fig. 6), the np differential cross 
section for forward angles decreases when going from p=0 to po and then increases 
for higher densities. The differential cross section at backward angles always decreases 
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with density. While the free np differential cross sections are highly anisotropic, the 
in-medium cross sections become more isotropic with increasing density. 

We have also calculated the in-medium cross sections using other potentials (Bonn B 
and C) and found very little difference as compared to Bonn A. Thus, there is fortu- 
nately little model dependence in our predictions. 

It is also interesting to compare the results obtained in this work with the parametrized 



NN cross sections proposed by Cugnon et al. p2|, ^ which are often used in trans- 



port models such as BUU and QMD. This comparison is included in Fig. 6. It is 
clearly seen that, while the Cugnon parametrization is almost isotropic, the micro- 
scopic results still have some anisotropy at all densities considered. The anisotropy 
in the present results decreases with increasing density. There is also density depen- 
dence in the microscopic differential cross section, while the Cugnon parametrization 
is density independent. 

We mention that Cugnon et al. have parametrized the free pp cross sections. 
This explains the almost-isotropy in their differential cross sections as well as the 

m 



lack of the density dependence. Note that in the work of Cugnon et al., no 
difference is made between proton and neutron; thus, the pp cross sections are also 
used for np scattering. There are, however, well-known differences between pp and np 
cross sections which in the more accurate microscopic calculations of the near future 
may be relevant. The difference between the (in-medium) np and pp cross sections 



is discussed in detail in Ref. |Q. In this section, we restrict our discussion of "NN 
cross sections" to np cross sections. 

The np differential cross section in free space can be well parametrized by the 
following simple expression: 

^{Eiab, 0) = ^^^^ -explhicosH + sm^- - 1.0)1 (106) 

dn^ ' 1.0 + 0.05(EfJ - 15.5) ^ 1 ^ ' 



with 
and 



0.0008(Ef„f - 4.625) for Eiai < lOOMeV 



h = 0.0006(36.65 - Ef^f ) for Eiab > lOOMeV. 

with Eiab in the units of MeV. 

It would be useful to parametrize the in-medium np differential cross section 
as well. However, the complicated dependence of the in-medium differential cross 
sections on angles, energy, and especially density makes this very difficult. Instead, 
we have prepared a data file, containing in-medium differential cross section as a 
function of angle for a number of densities and energies, from which the differential 
cross sections for all densities in the range 0-3po and all energies in the range 0-300 
MeV can be interpolated. This data file is available from one of the authors (R. M.) 
upon request. 

In addition to the in-medium NN differential cross sections which enter the trans- 
port models to determine the direction of the outgoing nucleons, the in-medium total 
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NN cross sections, a^N, are also of interest. They provide a criterion for whether a 
pair of nucleons will collide or not by comparing their closest distance to ^Jctnn/'^- 
We show in Fig. 7 the in-medium total cross sections as function of the incident en- 
ergy Eiab and density. It is seen that the in-medium total cross sections decrease 
substantially with increasing density, particularly for low incident energies. 

Figure 7 also includes a comparison with the total cross sections used by Cugnon 



et al. ^ (dotted line in Fig. 7). It is seen that at low energies and low den- 
sities, the Cugnon parametrization underestimates the microscopic results, while at 
higher energies and higher density it is the other way around. Note that the Cugnon 
parametrization is not density dependent, and, thus, predicts the same for all densi- 
ties. 

Finally, we propose a parametrization for the total np cross section as a function 
of the incident energy Eiab and density p 

10-1-0 nrr^zL pi-si n2 

anMat,p) = (31.5 + 0.092abs(20.2 - EUTY-^^y^^ (107) 

where Eiab and p are in the units of MeV and fm~^, respectively. 

The major conclusions of the present microscopic calculations are: 

• There is strong density dependence for the in-medium cross sections. With the 
increase of density, the cross sections decrease. This indicates that a proper 
treatment of the density-dependence of the in-medium NN cross sections is 
important. 

• Our microscopic predictions differ from the commonly used parametrizations of 



the differential and the total cross sections developed by Cugnon et al. [|82|, p2| . 
The Cugnon parametrization underestimates the anisotropy of the in-medium 
np differential cross sections. In the case of the total cross sections, the Cugnon 
parametrization either underestimates or overestimates the microscopic results, 
depending on energy and density. 



More details of this investigation can be found in Refs. |91|, |93| . 



7 Finite Nuclei 

Encouraged by the good results for nuclear matter, one may now try to describe 
finite nuclei starting from the free-space nucleon-nucleon interaction. A straightfor- 
ward way would be to solve the relativistic Brueckner-Hartree-Fock equations for 
finite systems. This is, however, an extremely difficult task [^, |98|. Therefore, it 



might be a reasonable next step to incorporate the DBHF results in a relativistic 
Hartree framework where the coupling constants are made density dependent so as 
to reproduce the nuclear matter results. This relativistic density dependent Hartree 
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(RDDH) approach is similar in spirit to the work by Negele in the nonrelativistic 
approach [^ . 

The working basis of the RDDH approach for finite nuclei is the relativistic Hartree 
Lagrangian ||29| (sigma-omega model Lagrangian of Walecka [0). Writing explicitly 
the density dependence of the coupling constants, we have 



-Crddh = - M - g^{p)a - g^li^^^) ^ 

+ 1 {d^af - \mla^ - \ {d,uj^ - d^uj,f + lm>J , (108) 



in conventional notation |Q. In Hartree approximation (mean field approximation) 
the nucleon self-energy in nuclear matter is given by 



^RDDH 



(p) = U,(p) + U,(p)7o. (109) 



Here, the scalar and the vector potentials are expressed in terms of the couphng 
constants go-(p) and guj{p) through 

U. p = —ps 

ml 

U.(p) = ^4^P., (110) 
ml 

where ps and p„ are the scalar and the vector densities, respectively: 

7 , . / d k M* 

^ ^ 
fcp 

(2^ ■ ^^^^^ 

The scalar potential is related to the effective mass by M* = M + Us(p). The 
connection of the RDDH approach to the DBHF theory is made through the nucleon 
self-energy in nuclear matter, Eq. (109). In fact, one can express the DBHF self-energy 
in nuclear matter in this form [^] using the symmetry requirements and redefinition 
of various terms through the use of the Dirac equation The density dependent 
coupling constants are then obtained through Eqs. (110) and (111), where Us(p) and 
Uv{p) are the results of the DBHF calculations of nuclear matter. 

We write the equations of motion for finite nuclei for completeness. The normal 
modes of the nucleon are calculated with the Dirac equation, 

(-za- V + /3M*(r) + F(r))^,(r) = EM^) (112) 

with M*(r) = M + g^{r)a{r) and V(r) = g^{r)uj^{r) + e^A°{r). 
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Figure 8: Vector and scalar potentials for relativistic density dependent Hartree (full 
lines) and relativistic Hartree only (dashed lines) as a function ofkp. 

The Klein-Gordon equations for a, uj^, and are 



where ps, Pv and pp are the scalar, vector and proton densities, respectively, which are 
obtained through the Dirac wavef unctions. We solve the coupled differential equations 
self-consistently. The center of mass corrections are applied in the same way as in 



In Fig. 8, the vector and the scalar potentials, U^, and U<j are depicted as a function 
of kp. At normal matter density, kp = 1.35 fm~^, U„ = 274.7 MeV and 11^ = - 
355.7 MeV. As a comparison, the vector and scalar mean-field potentials are shown 
using constant (density independent) couplings; i.e. the scalar and vector potential 
parameters are the ones of Potential A. This comparison clearly indicates that one 
needs density dependent coupling constants in order to reproduce the nuclear matter 
results within the relativistic Hartree framework. 

Now we fit the coupling constants go-(p) and ga;(p) to the self-energies and U„ 
with the choice that m^ and m^^ are the masses in Potential A as a standard choice; 
m^- = 550 MeV and m^^ = 782.6 MeV. As can be seen from Eq. (110), the relevant 
quantities are g^ and g^. At kp = 0.8, 1.1, 1.5 we obtain g^/47r = 12.3, 8.91, 6.23 
and g^/47r = 18.63, 13.48, 9.06, respectively. This shows that the coupling constants 




{-A + ml)uj\r) 
-AA°(r) 



9u^{r)p^{r) 



(113) 



Ref. 0. 
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at the surface are more than 40% bigger than in the interior. At nuclear matter 
density {kp = 1.35fm~l) the couphng constants are not modified as one can see from 
Fig. 8. But at smaller densities both coupling constants are growing. We note here 
that Us and U^, are, in principle, dependent not only on the density but also on the 
momentum of the nucleon. We neglect this weak momentum dependence in our study 
of finite nuclei. 

We calculate ^^O and ^"Ca as examples for finite nuclei within the RDDH ap- 
proach. The calculated results on the binding energies, single particle energies and 
charge radii (using the Bonn A potential) are tabulated in Table 9. For comparison, 
results of nonrelativistic Brueckner-Hartree-Fock calculations for ^^O with the same 



potential (Bonn A) are shown in the column N-BHF [^]. It is interesting to find that 
the RDDH results are very close to experiment. The improvements, as compared to 
N-BHF, are remarkable. The root mean square charge radius is almost perfect, while 
the binding energy is slightly smaller than experiment. If we compare relativistic 
and nonrelativistic results, we find that the radius of the relativistic calculation is 
larger. This is natural since, e. g., the lower component of the relativistic lp3/2 - 
wavefunction looks like a nonrelativistic ld3/2 - wavefuction. This shifts part of the 
density to the surface and leads to a larger radius in the relativistic case. This has 
also consequences for the binding energy per nucleon. Using the Bethe-Weizsaecker 
mass formula, we can estimate that surface and Coulomb effects lead to 2 to 3 MeV 
less repulsion for the relativistic calculation since the radius is larger. Although the 
volume effect is 1 to 2 MeV more repulsive in the relativistic case which has been 
determined in nuclear matter in Ref. the relativistic description of finite nuclei 
yields altogether more binding energy per nucleon. This shows that relativistic effects 
lead off the Coester band, which exists for finite nuclei, too. 

In Fig. 9, we show a comparison with experimental charge density distributions 
obtained from elastic electron scattering The RDDH results with the Bonn A 
potential (dash-dotted curve) compare fairly well with the experimental data. Look- 
ing more closely, the central density comes out to be higher than experiment and the 
density falls off slightly faster than experiment. This observation of the density and 
the slightly smaller binding discussed above, seems to be the reflection of the equa- 
tion of state of nuclear matter as shown in Fig. 4. In fact, when we take the DBHF 
results from the Bonn C potential, whose saturation density is almost perfect but the 
saturation energy is about 4 MeV above the experimental value, the charge density 
distributions are found to be very good as can be seen in Fig. 9 (dashed curve). In 
this case, however, the binding energies of flnite nuclei are found to be too small; i.e. 
E/A = - 5.9 MeV for ^^O and E/A = - 6.0 MeV for ^°Ca. 

Another issue is the density-dependence of meson masses; this is discussed in 
Ref. m. 
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Table 9: The binding energy per nucleon, root-mean-square charge radius and single 
particle energies are displayed for the relativistic density dependent Hartree(RDDH) 
approach, the nonrelativistic Brueckner-Hartree-Fock(N-BHF) method and experi- 
ment, respectively. The Bonn A potential is used. The upper part of the table is for 
^^O and the lower one for ^'^Ca. 



16Q 


RDDH 


N BHF 


Ex])eriiiicnt 


BE/ A [MeV] 


-7.5 


-5.95 


-7.98 


Yc [fm] 


2.66 


2.31 


2.70 ±0.05 


e(lsi/2) [MeV] 


-43.5 


-56.6 


-40 ±8 


e(lp3/2) [MeV] 


-21.8 


-25.7 


-18.4 


e(lpi/2) [MeV] 


-16.5 


-17.4 


-12.1 






RDDH 


N-BHF 


Experiment 


BE/ A [MeV] 


-8.0 


-8.29 


-8.5 


Tc [fm] 


3.36 


2.64 


3.5 


e(lsi/2) [MeV] 


-53.3 




-50 ±8 


e(lp3/2) [MeV] 


-36.0 




-34 ±5 


e(lpi;2) [MeV] 


32.5 




-31 ± 5 


e(ld5/2) [MeV] 


-19.3 


-30.2 


-14 ±2 


e(2si/2) [MeV] 


-14.3 


-24.5 


-10 ± 1 


e(ld3/2) [MeV] 


13.6 


-16.5 


-7± 1 



8 Summary and Outlook 

In this chapter we have presented the formalism for a relativistic approach to two- 
nuclcon scattering and nuclear structure. The latter is based upon a relativistic 
extension of Brueckner theory. The essential idea of this approach is to use the Dirac 
equation for the single particle motion. The nucleon self-energy in nuclear matter 
obtained in this framework consists of a large (attractive) scalar and (repulsive) vector 
field. The size of these potentials (several hundred MeV) motivates the use of the 
Dirac equation. 

Furthermore, we have applied relativistic meson-exchange potentials appropriate 
for this approach. The one-boson-exchange (OBE) model includes the six non-strange 
bosons with masses below 1 GeV/c^, 7i,rj,a,S,uj, and p. The potentials describe low 
energy NN scattering and the properties of the deuteron quantitatively. Thus, they 
are suitable for (parameter-free) microscopic nuclear structure calculations. Apart 
from the usual interaction Lagrangians for heavy mesons, these potentials apply the 
pseudovector (gradient) couphng for the ttNN (and rjNN) vertex. For a relativistic 
approach, it is necessary to apply this coupling for the pseudoscalar mesons, since the 
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Figure 9: Comparison of the densities in relativistic density dependent Hartree and 
experiment for ^^O and ^°Ca, respectively. The results with the Bonn A potential 
are shown by the dash-dotted curve, while the corresponding ones for Bonn C are 
depicted by the dashed one. 

commonly used pseudoscalar coupling leads to unrealistically attractive contributions. 
This has its origin in chiral symmetry. 

The consideration of the two-nucleon system is based upon the relativistic three- 
dimensional reduction of the Bethe-Salpeter equation introduced by Thompson. We 
derive and present the Thompson equation for an arbitrary frame. In two-nucleon 
scattering, this equation is applied in the two-nucleon cm. system, while in nuclear 
matter, the nuclear matter rest frame is used. Thus, wc calculate the nuclear matter 
G-matrix directly in the nuclear matter rest frame in which it is needed to obtain 
the energy of the many-body system. The advantage of this method is that, in the 
nuclear medium, the rather involved and tedious transformation between the two- 
nucleon cm. frame and the nuclear matter rest frame is avoided. 

Using the framework outlined, the properties of nuclear matter are calculated in 
the Dirac-Brueckner-Hartree-Fock (DBHF) approximation. The single particle scalar 
and vector fields as well as the single particle wave functions (Dirac spinors) are de- 
termined fully self-consistently. The strong attractive scalar field leads to a reduction 
of the nucleon mass in the medium, which increases with density. One consequence of 
this is a suppression of the (attractive) scalar-boson exchange. This strongly density- 
dependent effect improves nuclear saturation considerably such that the empirical 
saturation energy and density can be explained correctly. Our calculations include 
also the saturation mechanisms of conventional (non-relativistic) many-body theory 
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(i. e. Pauli and dispersion effects). The predictions for the compression modulus of 
nuclear matter as well as the Landau parameters are in satisfactory agreement with 
empirical information. 

In spite of the success of the present calculations, several critical questions can be 
raised. First, there will be contributions of higher order in the conventional hole-line 
expansion. For non-relativistic Brueckner theory. Day found an increase in the binding 
energy per nucleon of 5-7 MeV from the three- and four-hole line contributions. Note 
that in the work of Day the standard choice ('gap' choice) for the single-particle 
potential is used (this is true for all results displayed in Table 1 and Fig. 1).. However, 
in the calculations of our work the continuous choice for the single-particle potential, 
Eq. (98), is applied. In the two-hole line approximation, this choice leads to about 4 
MeV more binding energy per nucleon as compared to the gap choice. Investigations 
by the Liege group ||100|| suggest that a lowest order calculation with the continuous 
choice effectively includes the three-hole line contributions. In the light of this result, 
the contributions of higher order in the hole-line expansion, missing in our work, may 
be believed to be small. 

In a recent study, a special class of ring-diagrams has been summed up to infinite 
orders using our relativistic G-matrix. Only a minor change of the present DBHF 
result is found ||101|| which, however, further improves the saturation density. 

Secondly, one may question the OBE model. The weakest part of that model is 
the scalar isoscalar a boson. However, as mentioned, Dirac-Brueckner calculations 
have also been performed with more realistic meson models, in which the fictitious 
(J is replaced by an explicit description of the 27r-exchange contribution to the NN 
interaction. In those calculations, it is found that the relativistic effect comes out 
almost exactly the same as in the OBE model P3[ . This is not too surprising since it 



has been known for a long time that the 27r-exchange part of the nuclear potential is 
well approximated by the exchange of a single scalar isoscalar boson of intermediate 
mass. 

One may also criticize that the present approach is restricted to positive-energy 
nucleons. At a first glance, this may appear rather inconsistent in a relativistic 
framework. However, calculations done for two-nucleon scattering have shown that 
the contributions from virtual anti-nucleon intermediate states (pair-terms) are ex- 
tremely small (when the pseudovector coupling is used for the pion) Thus, 
the medium effects on this small contributions will be even smaller. This has been 
confirmed by the relativistic G-matrix calculations including antiparticle intermediate 
states by Amorim and Tjon [|102|| . 

Since higher orders in the hole-line expansion are attractive, while corrections from 
the Dirac sea (of anti-nucleons) have (as far as calculations exist) yielded additional 
repulsion |^, cancelations may occur between many-body contributions missing in 
the present approach. Thus, it is not unlikely that DBHF may ultimately turn out to 
be a reasonable approximation to this very complex (relativistic) many-body problem. 

We have also calculated cross sections for NN scattering in the nuclear medium. 
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We find strong density dependence for tlie in-medium cross sections. Tlie results from 
our microscopic derivation differ substantially from simple cross section parametriza- 
tions that are commonly used in heavy ion reaction calculations. Our in-medium 
cross sections together with our momentum-dependent mean field Q], which are 



both derived on an equal footing, can now be used as input for consistent transport 
model calculations of heavy ion reactions. 

Moreover, we developed a relativistic density-dependent Hartree approach for fi- 
nite nuclei, where the coupling constants of the relativistic Hartree-Lagrangian are 
made density-dependent and are obtained from the relativistic Brueckner-Hartree- 
Fock results of nuclear matter. The results on binding energies and root-mean-square 
radii of ^^O and ^°Ca agree very well with experiment. 

Further applications of our approach to the nucleon mean free path and proton- 



nucleus scattering are presented in Refs. ||103| , |104 . 

On a more fundamental level, one may raise the critical question how serious and 
reliable the meson model is. Notice that the relativistic effects obtained in Dirac 
approaches are intimately linked to the meson model for the nuclear force. In the 
meson-exchange model for the NN interaction, the isoscalar vector meson uj and the 
isoscalar a boson provide the largest contributions. In our many-body framework, 
they are responsible for the large effective scalar and vector fields obtained for the sin- 
gle particle potential (nucleon self-energy) in the many-body system. Assuming that 
quantum-chromodynamics (QCD) is the fundamental theory of strong interactions, 
one of the greatest challenges we have to face is the simple question: why do these 
meson models work so well — for two as well as many nucleons? Chiral symmetry 
may provide the key to the answer. 



This work was supported in part by the U. S. National Science Foundation under 
Grant No. PHY-9211607 and by the Deutsche Forschungsgemeinschaft (SFB 201). 
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